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Abstract 

A generalization of the non-Abelian version of the CP N ~ 1 models (also known 
as Grassmannian models) is presented. The generalization helps accommo- 
date a partial breaking of the non-Abelian gauge symmetry. Constituents of 
the composite gauge fields, in many cases, are naturally constrained to be- 
long to an anomaly free representation which in turn generates a composite 
scalar simulating Higgs mechanism to break the gauge symmetry dynami- 
cally for large N. Two cases are studied in detail: one based on the SU(2) 
gauge group and the other on SO(10). Breakings such as SU(2)— >U(1) or 
SO(10)— ►SU(5)xU(l) are found feasible. Properties of the composites fields 
and gauge boson masses are computed by doing a derivative expansion of the 
large N effective action. 

PACS nos: ll.15.Ex, 12.50.Fk 

I. INTRODUCTION 



Compositeness has been the way of nature. Some of the so called elementary particles 
of earlier times have turned out to be composites of more elementary ones. Compositeness 
is one way of getting to a simpler theory involving usually fewer fields or fewer parameters 



at the fundamental level. It is also expected to soften the ultraviolet behavior. Composite 
gauge fields could in addition provide us with an understanding of the gauge principle. They 
have received considerable attention in the literature [0-0] . The models that are relevant 
to our purpose are the so called CP^ -1 models and their non-Abelian generalizations 
called Grassmannian models(GMs) P|. Here the gauge fields arise as composites of bosonic 
constituents 0. These models have been studied in the large N limit and it is found that 
the non-Abelian symmetry is either completely broken or not broken at all. Ultimately one 
would like to construct phenomenological models along these lines, but the phase structure 
of these models is not very useful for that purpose. One needs a version in which the gauge 
symmetry is partially broken. 

Remarkably, as we will see in this article, there does exist a generalization of the GMs 
that allows for partial symmetry breaking. Some of the results of this paper have been briefly 
reported earlier ||. The primary agent of symmetry breaking turns out to be a scalar that 
too is composite. This composite Higgs scalar arises naturally as a solution of the modified 
constraint equation. In many of the cases, it belongs to the adjoint representation of the 
gauge group. The constituent fields in those cases belong to an anomaly free representation. 
One may recall here that that the agent of symmetry breaking in grand unified theories is 
usually a Higgs scalar in the adjoint representation. One may further recall that the fermions 
in a physical theory belong to an anomaly free representation, and we encounter the same 
feature here, though in the bosonic version. Two examples are studied to illustrate the 
approach: one based on the gauge group SU(2) and the other based on SO(10). The SU(2) 
example is the simplest and best suited to illustrate the approach. Here there exists a phase 
where SU(2) breaks to a U(l) subgroup. The case of SO(10) studied in some detail is more 
interesting from the physical point of view. The phase structure is richer with symmetry 
breaking to various subgroups such as SU(5) or SU(5)xU(l). 

We also compute the properties of the composite fields, the gauge bosons and the Higgs 
scalars, by doing a derivative expansion of the large N effective action. The expansions 
available in the literature do not serve our purpose as they are, to our knowledge, also 
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expansions in the Higgs scalar. We hence develop a suitable derivative expansion which 
we use to compute the kinetic terms and the mass terms for the composites in the various 
phases. 

To start with, in the next section, we review the known models relevant to our work. 
First we look at the CP N ~ X model that involves a U(l) gauge theory. Then we discuss its 
non-Abelian generalization, the Grassmannian model. It is based on the gauge group U(M) 



with the scalars in the fundamental representation. In section |TTT| , we introduce our model, 
a generalization of GM that is capable of accommodating other gauge groups with more 
general scalar representations. A suitable potential responsible for a rich phase structure is 
then introduced. In section |V|, we first illustrate our approach with a simple example based 
on the gauge group SU(2) and then look at the interesting but more complicated case of 
SO(10). Section [V| discusses the properties of the various composites. The global symmetry 
with its breaking patterns and the resulting Goldstone modes are discussed in section [Vl|. 
Section |V11| concludes with a discussion of the present approach. Derivative expansion of 
the effective potential useful in section |V] is carried out in appendix [A[ 



II. KNOWN MODELS 

The known models that are relevant to our purpose, the CP N_1 and the Grassmannian 
models, are briefly reviewed in this section. The simplest model is the one that induces a 
U(l) gauge theory, the CP N ~ l model. This is just a field theory involving N > 1 complex 
scalar fields, 

Z=(Z 1 ,Z 2 ,---,Z N ), (1) 

satisfying the constraint J2i = 1- The constraint sets their overall scale. A U(l) gauge 
invariance removes in addition an angular variable. Thus the model is in effect a field theory 
of N — 1 complex scalars. For convenience in writing the various results, we have above 
represented the Z^s collectively as a row vector Z. In this notation, the constraint can be 
rewritten as ZZ' = 1. The Lagrangian is 
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(5N 
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Here j3 is the inverse of a coupling constant. An overall multiplicative factor N is introduced 
for later convenience in the 1/N expansion. It is easy to verify the existence of a U(l) 
gauge invariance under which each of the Z^s transform with the same phase. In the vector 
notation, this is simply Z — * e Z where the phase 9 is space-time dependent. The constraint 
is clearly invariant under this symmetry. To see that it is a gauge symmetry, note first that 
the combination iZd^Z^ transforms as a U(l) gauge field, 
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zd^ -> izd^ + d^e. (3) 



The last term actually has a ZZ^ but that drops out due to the constraint. The U(l) gauge 
invariance will be more explicit, if we rewrite the Lagrangian by introducing an auxiliary 
field = iZd^Z^ as follows: 

L = (3N (d^Zd^ - 2iA li Zd IM tf + Al) = (3N [d^Z{D^] . (4) 

Here D^Z is the covariant derivative (<9 M — %A^)Z. In this form, the gauge symmetry is 
manifest. As shown in section |V|, the auxiliary field A^ that transforms as a gauge 
field becomes dynamical and hence a genuine gauge field after quantum corrections in the 
large iV approximation. It is a composite gauge field made of the Z fields. Thus the model 
under consideration can be viewed as an induced U(l) gauge theory or a theory of composite 
gauge fields. This model is a special case of the Grassmannian model; hence we study its 
phase structure below as a special case of GM. 

The Grassmannian model is a generalization of the CP^ -1 model that induces a non- 
Abelian gauge theory. We now have more fields, a set of them, represented collectively by 
a M x N matrix Z with the elements Z ai , a labeling the rows and % labeling the columns. 
The column index i is an internal index or a flavor index that is essentially carried over from 
our previous model. The new index a, the row index, is the gauge index associated with a 
non-Abelian symmetry which in the present case is U(M). All our results should reduce to 
those of the previous model for the case of M — 1. The constraint is now 



ZZ 1 * = I M , (5) 

where Im is an identity matrix of order M. The present Lagrangian is of the previous form 
(0), but now as such it will be an M x M matrix and hence needs an overall trace to make 
it a number, 



L = pNti 



(6) 



It is again easy to verify that there is a U(M) gauge invariance with respect to the index 
a. Under this gauge symmetry, Z transforms as a set of N fundamental representations. In 
the matrix notation, this transformation is simply Z — > UZ , U being an M x M space-time 
dependent unitary matrix representing the gauge transformation. This transformation does 
not affect the i index which labels N fundamental representations. The constraint respects 
this symmetry. The object = iZd^Z^ is in the adjoint representation of U(M) and 
transforms as a gauge field thanks to the constraint, 

A^UA^ + iUd^. (7) 

The role of the constraint here is to simplify the last term above from iUZZ^d^U^ to iUd^U'. 
The gauge symmetry becomes explicit when we rewrite the Lagrangian as in Eq. (f|) with 
an overall trace, 

L = pNti [D^D^] . (8) 

D^Z being the covariant derivative (<9 M — iA^)Z. As before, A^ appears as an auxiliary 
field but, as can be seen at large N, it becomes dynamical and hence a genuine gauge field 
after quantum corrections. It is a composite gauge field with the Z fields as constituents. 
The constraint and the U(M) gauge invariance have the effect of suppressing M 2 degrees of 
freedom. The theory is thus based effectively on M(N — M) scalars. Clearly, for it to be a 
sensible one, iV is required to exceed M. 

The constraint ZZ< = Im can be incorporated into the Lagrangian with the help of a 
Lagrangian multiplier E, a M x M matrix. The result is 



L = 0NtT iD^ZiD^Zy + EZZ ] - E 



(9) 



To understand symmetry breaking, and hence to identify the various phases, we need to 
obtain the effective potential. We will do this at large N. Because is not expected 
to pick up any expectation value, we will set it to zero. The classical contribution to the 
effective potential comes from the Lagrangian (Q) by dropping the derivative terms. Because 
1/N appears in the Lagrangian like the Planck's constant, the quantum corrections to this 
contribution is expected to be suppressed by a factor 1/N. But there are N fundamental 
representations contributing equally and this can offset the 1/N suppression. The result is 
that at large N the effective potential for the Z and £ fields obtained by integrating away 
the Z fluctuations carries a correction 
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Here and in the rest of this article, we suppress the dependence of the momentum integrals 
on a cutoff A. The total effective potential is thus 



V BS = (3Nti [EZZ^ -E)+Nj A^tr In (k 2 I M + s) . 



(11) 



To determine the various phases, we need to extremize this potential. The resulting 
saddle point equations (SPEs) are T.Z = obtained by varying Z^ and 
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(2tt) 4 m M + E 

coming from varying S. Now, let us look for solutions of the form 
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where I p and Im-p are t wo identity matrices of order p and M — p respectively. Note that 
EZZ^ is zero with this ansatz. Solutions for Z that satisfy = can easily be constructed. 
Eq. (|T2"D leads to two equations, one in the I p sector and the other in the Im-p sector, 
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d 4 k 1 
(2tt) 4 k 2 + a 
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(14) 



J (27r) 4 F 



For p = there is no J p sector and hence the first equation above would be absent. Similarly 
for p = M the second would be absent. Let us first show that p can not lie in between. We 
will do this by showing that the two equations can not be satisfied simultaneously. Firt note 
that cr should not be negative for the momentum integral involving it to be well-defined. 
Hence, from the first equation, we note that (5 has an upper limit f3 c given by 



But the second equation implies that (3 C is also a lower limit of (3. To see this, rewrite it 
as v 2 = 1 — Pc/P and note that v 2 can not be negative. It thus follows that, except at the 
critical point (5 — f3 c where the question is irrelevant, p is either zero or M. 

Thus we have two phases. For (5 > f3 c , we have the broken phase (p = 0) where ZZ^ has 
an expectation value. The solution for Z has an expectation value along all the "directions" 
in the fundamental representation. This breaks U (M) completely and all the gauge bosons 
are massive. For (3 < /3 C , we have the unbroken phase (p = M) where the gauge symmetry 
is unbroken and the gauge bosons are massless. (5 = f3 c is a critical point. In fact, (5 < f3 c is 
a critical line along which all the masses vanish. 

In other words, the gauge group is either completely broken or not broken at all. There 
are apparently no phases, at least at large N, where a partial breaking of the gauge group 
is possible. To obtain a richer phase structure, we invoke a generalization of these models 
and study them at large N in the following section. 



As discussed in the previous section, the Grassmannian model involves scalars in N fun- 
damental representations of the gauge group U(M). A natural extension is to construct 




(15) 



III. A GENERALIZATION 
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models for various gauge groups with the scalars transforming under different representa- 
tions. To our knowledge, they have not been studied in the literature. In this section, we 
explore an interesting class of these models that are nontrivial generalizations offering a rich 
phase structure. As usual, we are concerned with symmetry breaking at large N. 

A. Modifying the constraint 

One approach, a straightforward one, is to choose some other gauge group G in place 
of U(M) but to leave the constraint ZZ^ = Im unchanged. In other words, we take the Z 
fields to belong to an arbitrary representation R of dimension M and multiplicity N of some 
chosen gauge group G. We may still represent the Z fields in the form of an M x N matrix. 
The transformation matrix U is now in R acting on the matrix Z as before, Z — > UZ. The 
Lagrangian is still of the form we came across earlier in Eq. (|8|), but with the auxiliary gauge 
field now taking values in the Lie algebra of G. An expression for the auxiliary field and 
the form of the Lagrangian involving the Z fields alone can be easily derived. They are 
given respectively by the Eqs. fll?|) and (JTJ|) given below. However, symmetry breaking at 
large N remains the same. This is because the gauge fields are set to zero in our discussion 
of the phase structure. 

A more interesting generalization occurs when the constraint is modified as well. Again, 
we take the Z fields to be in any representation R of dimension M and multiplicity A of a 
gauge group G. We look for a Lagrangian that resembles (§). It is clearly gauge invariant 
with the auxiliary gauge field transforming as in ([?]). Note that the part of the Lagrangian 
quadratic in A^ = A^T a is proportional to 

A a Xtr(T a T b Ztf) = A a Xtr(T ab ZZ^), 

where T a 's are the generators of the gauge group G and T a b = (T a T{, + T&T a )/2. Earlier in 
the GM, the constraint ZZ^ = Im was responsible for rendering it quadratic in A^ alone. 
This resulted in a well defined expression for the auxiliary field A^ as a composite of the Z 
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fields. Now, more generally, we achieve the same goal by imposing the following constraint 
instead: 



tr(T ab ZZ^) = I5 ab . (16) 

Here / is the Dynkin index of the representation R defined by tr(T a T b ) = I5 ab . Note that this 
new constraint also respects the gauge symmetry. Using it, it is easy to obtain the following 
expression for the auxiliary field: 
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T a (Zd^ -d,ZZ^)}. (17) 



Later in section |V|, we observe that this field becomes dynamical at large N and is justly 
called a composite gauge field constructed out of the Z's. An expression for the Lagrangian 
in terms of the Z fields alone can now be obtained, 

L = pN\x(dpZd»tf) + j^N {tr [T a (Zd^ - d^ZZ^)] f . (18) 

This derivation should ensure gauge invariance and this is easily seen to be the case. Incor- 
porating the constraint into the Lagrangian using a Lagrange multiplier S = T, a bT a b leads 
to an expression that agrees with (^j. The symmetry breaking effects could be potentially 
different since the matrix £ is not an arbitrary one any more. 

First, we make sure that the constraint ([16]) is different from the earlier one ZZ^ = Im- 
ZZ^ = Im clearly is a solution of (|T6"|). Fortunately, there are cases where this is not 
the only solution. The new constraint is in some cases weaker than the earlier one. To 
see this, introduce a hermitian matrix W by ZZ* = Im + W and observe that the new 
constraint is equivalent to looking for a solution of tr(T a bW) = 0. Given a W that leads to a 
positive semidefinite ZZ\ Z is solvable generally as Z = (Im + WY^ 2 Z for some Z obeying 
ZqZq = Im- The earlier constraint corresponds to the trivial solution W = 0. That there 
exist cases where W is nontrivial can be seen as follows. Make the ansatz that W is in the 
Lie algebra itself, that is W — W a T a . Now, the constraint tr(T ab W) = simply states that 
the Adler-Bell-Jackiw anomaly associated with the representation R should vanish. That 
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there exist anomaly free representations is well known. A simple example is the doublet of 
the gauge group SU(2) that yields a triplet for W. We will use this example later to illustrate 
our approach. A more interesting example is the spinor representation 16 of the gauge group 
SO(10) that is well known to anomaly free. The ansatz gives a solution for W that is in 
the adjoint representation 45. This is interesting, for it is known that an adjoint scalar is a 
promising candidate to break a grand unified theory based on SO(10). Its appearance here 
is quite unexpected. 

In some instances, the above ansatz gives the most general solution. This is the case with 
both the examples mentioned above. This can be seen from representation theory. Regard 
W as belonging to the product representations 2x2 = 1 + 3 and 16 x 16 = 1 + 45 + 210 
in the SU(2) and SO(10) examples respectively. The component representations in these 
decompositions should separately obey the constraint. The singlet appearing in both the 
cases is ruled out easily. This is because the singlet component being proportional to identity 
gives ti(T ab W) oc 5 ab violating the constraint. In the case of SO(10), we need to exclude 
the 210 as well. It is easily observed that there exist 16 x 16 traceless hermitian matrices 
violating the constraint. This excludes 210 because such matrices can only have components 
along the 45 and 210, and the 45 alone can not violate the constraint. It is not to be deduced, 
however, that the ansatz always gives the most general solution. For instance, if one were to 
pick a sufficiently large representation of the gauge group for R, one will easily end up with 
more representations that remain unsuppressed in W. But, representation theory should 
still be applicable to solve for W in general. 

Let us call the models as type one models when W solves identically to zero and our new 
constraint reduces to the old one. They are closer to the Grassmannian models discussed 
earlier, or rather to their generalizations mentioned in the beginning of this section involving 
subgroups of U(M). The other models where W can be nontrivial is referred to as type two 
models. Models based on a reducible R are quite generally of type two. This is because the 
constraint does not determine some components of W, for instance those connecting different 
subrepresentations in R. Note, however, that a reducible R arising from an irreducible one 
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repeated say q times, though appears to be of type two, can be recast as of type one by 
combining q and N to an overall multiplicity qN in place of N. We prefer to view them 
as type one models with multiplicity qN. Type one models are thus necessarily based on 
irreducible i?'s. 

The effective potential at large N is again of the form ([□]) encountered earlier. The 
SPEs governing the phases are HZ = obtained by varying and 

/Jtr [T o6 (ZZt - /„)] + / jgjtr (l-j^) =0 (19) 

obtained by varying S. For type one models, the traces can be dropped and the equation 
becomes equivalent to the one we had earlier for GM (see Eq. (|P2])). Solutions can be sought 
in the same manner. The phase structure is governed by a critical (3 given by Eq. dl5[) . It 
is ensured that the ansatz for E one makes in solving these equations is consistent with its 
definition E = H a bT a b. This is because the final conclusion involves either the broken phase 
E = (for (3 > (3 C ) or the unbroken phase E = ol M (for (3 < (3 C ). This is equivalent to 
E a f> = or E a b = cr5 a b/C2(R), where C^-R) is the second Casimir invariant of the irreducible 
representation R defined by T a T a = Ci(R)Im, and is acceptable. 

For type two models, the above equation can still be reduced to that of GM, Eq. flT2D, 
but with a matrix W satisfying tr(T a bW) = on the r.h.s, 

The factor (3 on the r.h.s makes this equation agree with ZZ^ = Im + W at the level of 
expectation values. Again, the solutions E = (for (3 > (3 C ) and E = ol M (for (3 < (3 C ) will 
satisfy this equation, for W can clearly be chosen zero. As before, E = oIm is acceptable 
for irreducible i?'s as it follows from E ab = crS a b/C2{R). In the case of reducible i?'s, that 
is R — J2i Ri where each Ri is irreducible, the solution for (3 < (3 C is a bit more involved. 
Try again the ansatz E a 5 = <j5 a b- In this case E is not proportional to identity. Instead, it 
is a diagonal matrix taking values C2{Ri)cr along each representation R^. We look for a W 
matrix that is also diagonal, with values Wi along Ri. Note that the constraint tr(T a bW) = 
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requires the Wj's to satisfy J2ihwi = where the Zj's are the indices of the representations 
Ri'B. The SPE along Ri is 



This determines a. Individual equations simply determine the various Wj's. That this 
solution holds only for j3 < j3 c is easy to note. 

The solutions obtained so far are not a priori the most general ones. There could be more 
solutions. This does not, however, appear to be the case with the examples mentioned earlier 
based on the gauge groups SU(2) or SO(10) (see Eqs. (|28[) and ( |32"D and the discussions 
following them). Analogous situation occurs for the gauge group E e with its representation 
27. This is perhaps illustrative of a generic phenomenon or suggestive of the need to look at 
larger representations that might lead to more solutions. We do not wish to go into those 
details here, rather we find it more rewarding to consider the other possibility, that of adding 
a potential. 



The unexpected appearance of an adjoint scalar W apparently didn't help us in a par- 
tial breaking of the gauge group. The situation changes drastically when a potential is 
introduced leading to a rich phase structure. The adjoint scalar, that has not played any 
significant role so far, plays a major one in the presence of a potential. Note that there is 
no simple way to incorporate a potential in the canonical GM without spoiling the global 
symmetries or the constraint equation. But, interestingly, the generalized models of the 
previous section, governed by the gauge invariant Lagrangian of Eq. ([18]) constructed with 
scalars alone, do allow for potential terms. As we will see, the presence of a potential leads 




(21) 



The constraint on the lUj's gives 




(22) 



B. Adding a potential 
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to drastically different conclusions. These models with such phase structures are relevant in 
model building. 

Let us keep the potential quite general to begin with, f3NtrV(ZZ^), where V(-) is an 
ordinary function of its argument, a polynomial for instance. This is expected to be a 
nontrivial extension of type two models unlike the case of type one models in which the 
constraint ZZ^ = Im reduces this to the addition of a constant. It is convenient to introduce 
a composite field variable X for ZZ^ and write the potential as PNtrV(X). The requirement 
X = ZZ^ can be incorporated with the help of a Lagrange multiplier Y, adding a term 
/3Ntr(Y ZZ^ — YX) to the potential. As before, constraint fll6|) can be accommodated with 
the help of a Lagrange multiplier E = E a ftT a fc. Its effect is, as we know, to add a term 
(3Ntr(EZZ^ — E) to the potential. After translating Y to Y — E for convenience, the total 
Lagrangian looks like 

L = (3Ntr [D^Z(D^Zy + V(X) + YZZ X — YX + EX — s] . (23) 

The large N effective potential is now computable, 

V cS = PNti [V{X) + YZZ^ - YX + EX - e] + N J ^i tr ln ( fc2/ M + Y ) ■ (24) 

The SPEs are obtained by extremizing this potential. Varying X determines Y to be S + 
V'(X) where a prime denotes differentiation with respect to the argument. Varying E gives 
tr[T a t>(X — Im)] — 0- As before, one may look for a solution of this in the form X = Im + W 
where W satisfies tr(T a f ) W / ) = 0. Varying Y and using these solutions yields 

^ ZZ '- J ")+/(gl WM + S+ V, (/M + H ^ ty - (25) 
This is to be supplemented with YZ = [E + V'(Im + W)] Z = obtained by varying Zl 



This system of equations resembles the ones obtained earlier (see Eq. (^)), with E replaced 
by E + V'(Im + W). The presence of V'(Im + W), however, is suggestive of a different phase 
structure. 

For type one models W = 0, and V'(Im + W) just adds a constant to E. This can 
be absorbed into E because these models, being based on an irreducible R, allow for the 
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addition of term proportional to identity to S. As expected in the beginning of this section, 
this is a trivial extension. However, this is not the case for type two models and we expect 
a rich phase structure. 

An example we will use later is a potential of sixth degree in Z and Z\ that leads to 
V'(Im + W) = alu + bW + cW 2 for some constants a, b and c. For a model based on an 
irreducible R, the term clIm can be absorbed into £ as we have already noted. When W is 
in the adjoint representation, the term cW 2 can also be absorbed into S. As a result, the 
SPEs to be solved are 

and [S + bW] Z = 0. These equations are difficult to handle analytically and we present our 
numerical results for SU(2) and SO(10) below. We find that they do have solutions for a 
range of parameters when {3 < (3 C and b < 0. 

Given all such solutions for (3 < /3 C , the next step is to determine those that are preferred 
energetically. The GM solutions of section [TJ leave the gauge group unbroken for (3 < (3 C 
whereas those found here break it at least partially. Which one is preferred is of course 
determined by the effective potential. In other words, one needs to compute V e Q for all the 
solutions and pick the one (or more) that has the the lowest value. We do this along a 
chosen path in the parameter space of (3 and b that crosses all the phases. We find that 
some of the new solutions end up always having the lowest potential. In other words, for a 
range of parameters, a partial breaking of the gauge group is preferred over the unbroken 
case. Details are presented in the next section. The following expression for the effective 
potential at a saddle point (SP) is used to this end: 

Kff(SP) = -PNtr (S + bW 2 /2) +N J -0^ti In (k 2 I M + E + bW) . (27) 

We are, as always, concerned with an adjoint W for an irreducible R. In the above expression, 
terms al and cW 2 have been absorbed into £ for convenience. 
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IV. TWO EXAMPLES: SU(2) AND SO(10) 



A. The case of SU(2) 

The example based on SU(2) with doublet Z's is the simplest and the most convenient 
one to illustrate the ideas presented above. The matrix £ = £ a bT a 6 is now proportional to 
identity, hence chosen to be al 2 - The W scalar, being a triplet, is taken to be along the 
03 direction, that is W = wa^. First we consider the case when the Z fields develop no 
expectation value. The resulting SPEs in the presence of a potential are 

d A k 1 



r ark i 
~ P + J (2txY k 2 + a + bw ~ ^' 



,p + / -^4r- —— = -f3w. (28) 

H J 2n) 4 k 2 + a-bw H K J 



(2?r) 4 k 2 + a + bw 
drk 1 
(2tt) 4 k 2 + a-bw 

Note that there is no solution to these equations when 6 = 0, that is in the absence of 
a potential, other than the one discussed earlier where W = and E is proportional to 
identity. In the presence of a sixth degree potential given by V (Im + W) = alu + bW + cW 2 , 
these equations do have solutions for a range of parameters. This can be seen by treating 
x — (a + bw)/A 2 and y — (a — bw)/ 'A 2 as independent variables to determine (3 and w from 
the above two equations. Given x and y and knowing w, one obtains b from x — y = 2bw/A 2 . 
The resulting equations are 

Pc b x — y 

'^A 2 = x ln(l + TJx) -y ln(l + TJyj' ^ ' 

For the momentum integrals to remain well defined, x and y should be positive (or zero). 
The region of the parameter space of (3 and b is obtained by letting x and y vary from 
zero to infinity. This falls in between the curves (a) and (b) shown in Fig. [I]. There are 
two solutions for a given (3 and b in this regime, but they are related to each other by an 
interchange of x and y and should be treated as one. Note that all these solutions yield 
(3 < (3 C and b < 0. Curve (a) has x = y and is the critical line. In fact, the region below 
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curve (a) is a critical surface where all the masses vanish. Curve (b) has one of x, y zero. 
Symmetry breaking involved here is from SU(2) to U(l). 
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FIG. 1. The phase diagram in the case of SU(2) obtained solving Eqs. (29) and (31). The 
effective potential of Fig. 2 is computed along the dotted line. Parameters in the theory are (5 and 
b, and A is the momentum cutoff. Details are given in the text. 

There exist new solutions for a nonzero Z as well. Giving an expectation value diag(0, v 2 ) 
for ZZ\ the second equation in (p8|) gets replaced by 



where we have set y = a — bw to zero to satisfy (S + bW)Z = 0. Treating x = (a + bw = 
2bw)/A 2 and y' = v 2 f3/f3 c as independent variables to determine the others, 



one notes the presence of solutions in a parameter range for positive x and y' . Here as well, 
we require (3 < (3 C and a negative b. The parameter range is the one above curve (b) in Fig. 
|l|. SU(2) symmetry is now completely broken. Giving an expectation value diag(t> 2 ,0) for 



(3v 2 - (3 + (3 C = -f3w, 



(30) 




f3 c b x 



(31) 



(3 A 2 x ln(l + l/x) +y 
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ZZ^ is equivalent to this case and leads to no new solutions. A nonzero ZZ^ of the form 
diag(vj, v%) requires v\ = v\ and coincides with the completely broken case of GM discussed 
in section [D]. 

These are all the solutions. Now, consider all those for (3 < j3 c . The corresponding one of 
GM leaves the gauge group unbroken whereas those of this section break it at least partially. 
As discussed earlier, the effective potential needs to be examined to determine the preferred 
solution. We have chosen a path suitably fixing y crossing all the curves, shown dotted in 
Fig. |l|. Fig. ^| is a plot of the effective potential. The upper curve is for the GM solutions 
and the lower one is for the new solutions. Note that the lower sheet ends up always having 
the lowest potential. In other words, for /3 small, a partial breaking of the gauge group is 
preferred over the unbroken case. 




FIG. 2. A plot of the effective potential (with its zero appropriately chosen) versus (3 for a 
path shown dotted in Fig. 1 crossing all the curves. The crossings are denoted by (a) and (b). 
The upper curve corresponds to the unbroken case and the lower one corresponds to symmetry 
breaking as discussed in the text. 
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B. The case of SO(10) 



We now come to the second example, the gauge group SO(10) with the Z's in the 
representation 16. Here the possibilities for symmetry breaking are too many. We do 
not hope to address all of them, rather simply pick one possibility: SO(10) breaking to 
SU(5) or to its maximal subgroup SU(5)x U(l). Under SU(5), 16 of SO(10) decomposes 
to 10(1) + 5(— 3) + 1(5) where the U(l) charges are given in paranthesis. The unbroken 
symmetry for Z = corresponds to those generators that commute with the ansatz for E. 
To have symmetry breaking to SU(5)xU(l), we hence choose Z = and take S ab = a8 a b 
along the SU(5) directions, p along the U(l), zero otherwise. Note that a or b index runs over 
the adjoint representation 45 of SO(10) that under SU(5) decomposes to 24 + 10 + 10+1. Our 
ansatz for S afe corresponds to having it nonzero for (a, b) along (24, 24) and (1, 1). One could 
have it nonzero along (10, 10) and (10, 10) as well, but it turns out that this can be absorbed 
into a and p. This means that the S matrix is diagonal with values C2(10)<7 + p, C2(5)a + 9p 
and 25p along the representations 10, 5 and 1 respectively. With C 2 (10)/C 2 (5) = 3/2 and a 
suitable scaling of a, we may take them to 3a + p, 2a + 9p and 25p. The W matrix is taken 
to be along the U(l) direction; in other words, it is diagonal with values w, —3w and 5w. 
It is now straightforward to write down the SPEs in the presence of a potential, 



Consider first the case of no potential, that is b = 0. There are no new solutions. This is 
because the S eigenvalues 3a + p, 2a + 9p and 25p are either in the ascending order or in 
the descending order (or equal) and hence one cannot obtain the alternating signs on the 
r.h.s above. Allowing for nonvanishing Z does not improve the situation. Note that ZZ^ 
and £ should have nonnegative eigenvalues and T,ZZ^ = requires at least 3a + p or 25p 
to vanish to allow for a nonzero eigenvalue of ZZ^ . 




(32) 
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Solutions exist in the presence of a potential for a range of parameters. To see this, as 
in the SU(2) case, treat x = (3cr + p + bw)/A 2 and y = (2a + 9p — 3bw)/ A 2 as independent 
variables to determine (3 and w from the first two equations and z = (25p + 5bw) / A 2 from 
the last one. Given x and y and knowing z and w, one obtains b from 2x — 3y + z = 16bw/A 2 . 
Thus Eq. (^) yields 

3 3 ( 1\ 1 / l\ 

1 -^ = r ln ( 1 + x) + 4^H 1 + yj 
= (2x-3y + z)/4 

/3A 2 x ln(l + l/x) ln(l + l/y)' 1 j 

where z is a solution of 

z In (l + = 2x In (^1 + ~) - y In ^1 + ^ . (34) 

Note that x, y and z are required to remain positive (or zero) to keep the momentum 
integrals well defined. Solutions exist in a certain domain of x and y giving rise to a range 
for the parameters for 3 < (3 C and b < 0. The results of our numerical investigation is 
presented in Fig. [| There are in fact two solutions for a given 3 and b in the region between 
the curves (a) and (b), and one solution between the curves (b) and (c). In other words, 
one of the solutions extends from curve (a) to curve (b) while the other from curve (a) to 
curve (c). Curve (a) has x = y = z and is the critical line. Here too, the region below curve 
(a) is a critical surface wherein all the masses vanish. Curve (b) has z = and curve (c) 
has y = 0. The symmetry breaking involved here is from SO(10) to SU(5)xU(l) as noted 
before. 
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FIG. 3. The phase diagram in the case of SO(10) obtained solving Eqs. (33) and (36). The 
effective potential of Fig. 4 is computed along the dotted line. Parameters in the theory are (3 and 
b, and A is the momentum cutoff. Details are given in the text. 

Solutions exist for nonzero Z. Consider giving an expectation value v 2 for ZZ^ along 
the singlet in the decomposition 16 = 10 + 5 + 1. In this case, the third equation in (|32|) 
gets replaced by 



(3v 2 - (3 + p c = h(3w, 



(35) 



where we have set z = (25p + 5bw)/A 2 to zero to satisfy (S + bW)Z = 0. Again, treat x and 
y as independent variables and determine the others to obtain new solutions for positive v 2 
in a range of parameters: 



where we require 



& b = (2x- 3y)/4 

13 K 2 x ln(l + l/x) -y In (1 + !/?/)' 



^-v 2 = -2x \n(l + -) +y In f 1 + - J > 0. 



(36) 



(37) 



Pc V xj \ y / 

Here as well, one has (3 < (3 C and b < 0. The region of the parameter space covered by 
these solutions (one solution for a given (3 and b) is that in between the dashed curves (b) 
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and (d) of Fig. ^. Curve (d) has x = 0. The surviving symmetry here is SU(5) because 
a nonvanishing Z along the singlet in the decomposition 16 = 10 + 5 + 1 breaks the U(l) 
subgroup of SU(5)xU(l) as well. There are other possibilities. Giving an expectation value 
for ZZ^ along 5 (instead of the singlet) also leads to a solution which falls above the curve 
(c); solutions are also noted to exist for a nonzero ZZ^ along 10 and 1 (extending above 
curve (d)), or 10 and 5 (extending beyond that of 5). All of these, however, break the gauge 
group completely. 

What we have in the end is a two sheeted cover of the parameter space above the critical 
curve (a) in Fig. ^. One of them (call it the upper sheet) is through the solid curves while 
the other one (call it the lower sheet) is through the dashed curves. They meet along curve 
(a). There is of course one more sheet (call it the top sheet) for the solutions of our earlier 
case of the unbroken gauge group covering all of the parameter space for ft < ft c . This too 
meets the other two sheets along curve (a). For every point on any one of the sheets, there 
is a solution. 

As we have noted earlier, there could be more solutions. For instance, there is the possi- 
bility that a solution breaking SO(10) to SU(3)xSU(2)xU(l) (perhaps, with an additional 
U(l)) exists. The number of variables and the number of equations at least match, each 
being six, but the number of equations makes analysis complicated. There are more possi- 
bilities such as SO(10)— ►SU(4), SU(4) xU(l), etc. Each case has to be handled separately; 
a general treatment has eluded us. 

As before, which solution is preferred is determined by the effective potential. For this, 
one needs to compute V e s for all the solutions and pick the one (or more) that has the the 
lowest value. In the present case, this is not an easy task given the number of possibilities 
involved. Hence we will be content with doing this numerically for the solutions found above. 
We have chosen a path suitably fixing y in the lower sheet and z in the upper sheet crossing 
all the curves, shown dotted in Fig. |^. Fig. |] is a plot of this effective potential for the 
three sheets involved. The uppermost curve is for the top sheet, the middle one is for the 
upper sheet and the lowermost one is for the lower sheet. Note that the lower sheet ends up 
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always having the lowest potential. In other words, for (3 < f3 c but not close to it, a partial 
breaking of the gauge group is preferred over the unbroken case. 




FIG. 4. A plot of the effective potential (with its zero appropriately chosen) versus [3 for a path 
shown dotted in Fig. 3 crossing all the curves. The crossings are denoted by (a), (b), (c) and (d). 
The uppermost curve corresponds to the unbroken case and the lower two correspond to symmetry 
breaking as discussed in the text. 



V. PROPERTIES OF THE COMPOSITES 

Here we study the various composites encountered earlier, the gauge bosons and com- 
posite Higgs particles. Their properties can be read off from the effective action. Like the 
effective potential, the effective action at large N is obtained by integrating away the Z 
fields. It involves the original action plus a correction of the form 

S eff = -N Tr In [-D 2 + s) , (38) 

where 'Tr' stands for a complete trace, that is, a trace over the internal indices and an integral 
over the space-time coordinates. D is as before the covariant derivative = <9 M — iA^. For 
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simplicity, we have absorbed the bW term into E. To identify the kinetic terms for the 
various fields, a derivative expansion is needed. This is carried out in Appendix [A]. The 
result is that the kinetic terms for the S fields arise from 

f/ A /(0rf e -"7„'^^- <, -" , ^-" s < 39 > 

where D^O is [_D M , O] for any O. At the next order in the expansion, one obtains an 
analogous result for the gauge fields: 

These general results are valid for any gauge group. 

The two phases of GM, the completely broken and unbroken ones discussed in section 
U can be handled in this generic setup. In the unbroken phase where S = al, Z = the 
kinetic terms for the gauge fields simplifies to 

1 



2<? 2 (a/A 2 ) 
where g 2 is the coupling constant: 

1 N r d A k 



d 4 x tr F 2 (41) 



g 2 (a/A 2 ) 6 7 (2tt) 4 (P + a) 2 
iV 

ln(l + A7a 



96tt j 



1 + a/A 2 



(42) 



This also follows from the well-known contribution to the running gauge coupling constant 
at one loop order in the presence of N fundamental scalars. In the broken phase where 
£ = 0, ZZ^ = (1 — P c /P)Im, the above computation of the induced kinetic terms for the 
gauge fields suffers from an infrared divergence. Introducing an infrared cutoff \x for k 2 /A 2 , 
we obtain a result that coincides with the above one with a/A 2 replaced by fi. The mass 
terms arise from the kinetic terms for the Z fields that reads 

(3Ntr (AlZtf) = Nifi - &)tr (A 2 ) . (43) 

The mass squared for is then N{(3 — /3 c )g 2 (n). As expected, it is positive for (3 > f3 c and 
vanishes at the critical point (3 = j3 c . 
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We now turn our attention to type two models, in particular to the models based on 
SU(2) and SO(10). The mass squared results we obtain are all expressible in units of A 2 . 
All our results for the kinetic terms and mass terms involve the expectation values of the 
scalar field £ through x and y (and in the case of SO(10) z as well) parameters. In the 
phase where ZZ^ ^ 0, the expectation value of ZZ^ also appears in the expressions through 
v 2 . As we have seen in the previous section, all these parameters are expressible in terms 
the basic ones, (3 / j3 c and b/A 2 by solving a set of equations. When solving the equations 
however, we found it convenient to treat x and y as independent variables to determine 
P/Pc an d b/A 2 . This helped us to discover multiple solutions; but once we have chosen the 
energetically preferred solution, the relation between the x, y parameters and the basic ones 
is one to one and is hence invertible. 



A. The case of SU(2) 

First we consider SU(2) and compute the induced kinetic terms for the gauge fields in 
the U(l) phase where E = al 2 + bwa 3 and Z = 0. Writing 

Filv = 72 F/>+ + + \ F ^ azl (44) 

where a± = (oi ± «er 2 )/2, we obtain the kinetic terms for the gauge bosons of the unbroken 
U(l) and for those of the broken generators from Eq. (|4T)|). We have 

tr V^F^e-* 113 = e~ to j^cosh^) (i^) 2 + cosh[(t - 2t 1 )bw)F+F^ • (45) 

The kinetic term for the unbroken U(l) turns out to be 

1/1 I s 



d A x(Fl) , (46) 



8 V^ 2 (^) g 2 (y) / 

where g 2 is given in ( [42]) and as in section fV| x = (a + bw)/A 2 , y = (er — bw)/ A 2 . For the 



broken generators, the result is 

1 



2G 2 (x,y) 



d 4 xF^ u F~ u (47) 
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where 



N 



G 2 (x,y) 2(bw) 2 
The integral over k 2 can be done, 
1 N 



d A k 



k 2 + a, k 2 
In 



cr + bw 



2bw k 2 + a — bw 



(48) 



G 2 (x lV ) 87r 2 (x-y) 2 
where the functions I and J are 



x + y\ I(x) - I{y) l Jjx) - J(y) 

3 



x-y 



x-y 



(49) 



I(x) = ln(l + x) - x 2 \n(l + l/x) + x, 

J(x) = ln(l + x) + x 3 ln(l + l/x) + ar/2 - x 2 . 



(50) 



G 2 (x,y) is positive for x, y positive which is the region of interest. It is a monotonically 
increasing function of x + y in this region. It tends to g 2 (x) as we approach the critical 
line x = y. This is as it should be, since the SU(2) gauge symmetry is unbroken along the 
critical line and the kinetic terms for the broken and unbroken generators should add up to 
form SU(2) invariant kinetic terms. 

The gauge bosons of the broken generators receive mass terms. This arises from Eq. 
(BDI). Taking S to be space-time independent and writing 



(51) 



one first finds 



tr Xyr ( *~* l)E IV-* lS = 4e~* CT sinh[(t - ti)6^]sinh(ti6w)A+A M . 



(52) 



This yields for the mass terms 



where 



V 2 (x,y) = N 



d A k 



V\x,y) d 4 xA+A~ 



k 2 + a 



(53) 



1 , k 2 + a + bw 
In 



(k 2 + a) 2 — (bw) 2 2bw k 2 + a — bw 



(54) 
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Here too, the integral over k 2 can be done, but it is instructive to rewrite the result in terms 
of G 2 (x,y). Note that this integral can be obtained by differentiating the one in Eq. ([48]) 
with respect to a. This results in 

f 



(55) 



G 2 (x,y)_ 

It then follows that V 2 (x, y) is positive for x, y positive, the region we are interested in. As 
we approach the critical line x = y, 

t^)- (56) 

The mass-squared for is given by 

M 2 = G 2 (x,y)V 2 (x,y) 

^A 2 (x - y) 2 -^\ng 2 (x) as x -> y. (57) 

It is thus expressible in terms of the running coupling constant close to the critical line and 
vanishes along the critical line as expected. 

For the completely broken phase S = diag(a;, 0)A 2 and ZZ^ = diag(0, v 2 ), we still have 
the above results but with y = 0. The induced kinetic term for the U(l) field now suffers 
from an infrared divergence and suggests introducing an infrared cutoff /i for y. The mass 
squared result obtained above, though relevant with y = 0, now receives an additional 
contribution from the kinetic terms for the Z fields: 

j3Ntr (AlZtf) = l -PNv 2 A+ t A- + -^Nv 2 (4) \ (58) 

This makes the U(l) field A^ massive with a mass squared ~ /3Nv 2 g 2 (n)/2. The additional 
contribution to mass squared for A^ is j3Nv 2 G 2 (x, 0)/2. 

B. The case of SO (10) 



The computations for SO(10) are along the same lines. Consider the SU(5)xU(l) phase. 
Here £ = diag(x, y, z)A 2 along 10, 5 and 1, where as in section [TV] x = (3cr + p + bw)/A 2 , y = 
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(2(7 + 9p — 3bw)/A 2 , z = (25p + 5bw)/A 2 . One then easily computes the kinetic terms for 
the SU(5) generators: 



3 1 

+ 



(59) 



_g 2 (x) g 2 (y)_ 

where A runs over the 24 generators. The function g 2 has been defined earlier in Eq. (|42|). 
The part involving g 2 (x) arises from the 10, while the one involving g 2 (y) from the 5. Factor 
3 is a consequence of the fact that tr(T a T b ) for the 24 T's of SU(5) is 3 times in the 10 as 
in the 5. For the broken generators along the 10 and 10 appearing in the decomposition 
45 = 24 + 10 + TO + 1, we have 

3 1 



d x 



pa 



(60) 



_G 2 (x,y) G 2 (x,z) 

where a runs over the 10 generators and G 2 is the same function defined earlier in Eq. ( [4"9] 
The kinetic term for the U(l) field is 



1 

16 



2 9 5 



(61) 



_g 2 (x) g 2 (y) g 2 (z)_ 

where the superscript 45 denotes the U(l) direction. Note that as we approach the critical 
line x = y = z, gauge symmetry breaking disappears, G 2 — > g 2 and Eqs. (|5^), ( p0| ) and (|6TJ) 
add up to form SO(10) invariant kinetic terms as expected. 

It is straightforward to compute the mass terms for the 10 and 10 gauge bosons. First 
we note that 



tr V^'^V^e' 1 ^ = 4{3e"' (:E+s/)/2 sinh[(t - h){x - yV^sinh^a; - y)/2] + 

e-* (x,+2)/2 sinh[(t - h){x - z)/2]smh[t 1 (x - z)/2\) 



A" 
A* 



(62) 



where a A 2 has been absorbed into the t's in the r.h.s. Comparing this with Eq. (p2|) , we 
get for the mass terms 



3V 2 (x,y) + V 2 (x,z) 



A" 



(63) 



where V is the same function defined earlier in Eq. (p5|). This gives the mass-squared 
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M 



+ 



_G 2 (x,y) G*(x,z) 
for the 10's. As x,y and z tend to be the same, 

1 



3V 2 (x,y) + V 2 (x,z) 



(64) 



M 



8 



3(x — y) + (x — z) — lng (x). 

(JjJb 



(65) 



Here too, one obtains an expression in terms of the running coupling constant close to the 
critical line. It vanishes along the critical line as expected. 

For the SU(5) phase E = diag(a;, y, 0)A 2 and ZZ^ = diag(0, 0, v 2 ), the above results are 
still relevant but with z — 0. As in the case of SU(2), the induced kinetic term for the U(l) 
field suffers from an infrared divergence and suggests introducing an infrared cutoff ji for z. 
The mass squared result obtained above, though relevant with z — 0, receives an additional 
contribution from the kinetic terms for the Z fields: 



(3Ntr (A\ZZ 



\i5Nv 2 



A a 



(66) 



This makes the U(l) field massive with a mass squared f3Nv 2 g 2 (fx)/2. The additional 
contribution to the mass squared for the 10's is 

-l 



-/3Nv 2 



G*(x,y) G 2 (x,0) 



(67) 



VI. GLOBAL SYMMETRY AND THE GOLDSTONE MODES 

All the models we discussed have a global U(iV) symmetry. In addition to the gauge 
symmetry, this global symmetry could also suffer breakdown. It remains to investigate this 
breaking and the resulting Goldstone bosons and other massless particles if any. 

In the Grassmannian model, of the two phases, the unbroken phase retains this global 
symmetry. In this phase, only £ gets expectation value; but S is a singlet under the 
global symmetry. E expectation value however gives mass to all the Z scalars. The model 
has no massless particles in this phase. In the broken phase, Z gets an expectation value 
breaking the global symmetry in addition to the gauge symmetry. E expectation value is 
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now zero and all of the 2NM real components of Z are hence massless. Some of these are 
the would-be Goldstone bosons, eaten away by the broken gauge generators. Because the 
gauge symmetry U(M) is completely broken down, they are M 2 in number. It turns out 
that all those remaining, 2NM — M 2 in number, are the Goldstone bosons associated with 
the broken generators of the global symmetry. There are no unaccounted massless particles. 
To see this, choose the Z expectation value to be of the form 



where Im is an identity matrix of order M and Om,n-m is a zero matrix of order M x (N—M) . 
This expectation value breaks the global symmetry from U(iV) down to U(iV — M). The 
number of broken global generators are now easily computed; they are 2NM—M 2 in number. 
Along with the would-be Goldstone bosons, they account for all the 2NM massless particles. 

Coming to our type two models, the global symmetry is broken down when Z gets an 
expectation value. In the SU(2)— >U(1) phase or the SO(10)— ►SU(5)xU(l) phase of our 
examples, the global symmetry remains unbroken. E, a singlet under the global symmetry, 
also picks up an expectation value in this phase making all the Z scalars massive. There are 
no massless states. In the other interesting phase of our examples, Z picks up an expectation 
value along some direction, a singlet of SU(5) in the case of SO(10). £ expectation value 
in that direction is forced to zero. This will introduce 2N real massless states of Z. The Z 
expectation value can be arranged to be of the form 



where v is a column vector pointing in the singlet direction. This implies that the global 
symmetry is broken down from U(iV) to U(iV — 1). There are 2N — 1 Goldstone bosons 
associated with this breaking. The remaining one massless state of Z is a would-be Goldstone 
boson eaten away by the broken gauge generator. This is consistent with the fact that Z 
expectation value of the above type breaks one additional gauge generator. Again, there are 
no unaccounted massless states. 



Z oc (Im, Om,n-m) , 



(68) 



Z oc (v,O m ,n-i) , 



(69) 
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VII. DISCUSSIONS AND CONCLUSIONS 



We have presented an approach to composite gauge bosons that allows for partially 
broken gauge symmetries. It is a generalization of the well-known Grassmannian models 
that otherwise allow for either unbroken or completely broken gauge symmetries. In our 
approach, it is also possible to incorporate interesting potential terms leading to a rich 
phase structure. Even the simplest model based on SU(2) is not amenable to analytical 
handling of its phases; numerical investigation is called for. For models that are physically 
interesting in connection with unified theories, even a numerical analysis of all the phases is 
a challenging endeavor. 

We have illustrated our approach with an SU(2) example and analyzed in some detail 
an SO(10) example that could be of interest to unified models. What is remarkable of 
this exercise is that a set of equations governed by only two parameters gives rise to a 
rich set of solutions with interesting symmetry breaking patterns. There exist regions of 
the parameter space where SU(2) breaks down to U(l). In the case of SO(10), symmetry 
breaking to SU(5) or to SU(5)xU(l) or perhaps to some other subgroups is possible. These 
examples help realize our goal of constructing an induced gauge theory with composite gauge 
bosons having partial symmetry breaking. 

We have computed the properties of the composite fields, the gauge bosons and the 
Higgs scalars, by doing a derivative expansion of the large N effective action. Because we 
need an expansion that does not perturb the Higgs field, we cannot utilize the canonical 
expansions available in the literature. We have developed a suitable derivative expansion 
in the Appendix and have used it to compute the kinetic terms and the mass terms for the 
composites in the various phases. 

We have not addressed the issue of renormalizability of Grassmannian models or our 
generalized ones. It is interesting to note that the theory at large N exhibits a critical point 
which extends to a critical line in the presence of a potential. It is known that the critical 
points or lines can, and in many cases do, soften the ultraviolet behavior. This softening is 
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probably not sufficient enough to help renormalize the theory in four dimensions and one 
may have to include other relevant operators in the Lagrangian. In this connection, we note 
that certain four dimensional Grassmannian models of composite gauge fields have been 
studied on the lattice and shown to be renormalizable Their phase structures and their 
relation to continuum theories remain unexplored. 

This work is supported in part by U. S. Department of Energy, Grant No. DEFG-ER91- 
40672. 
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APPENDIX A: DERIVATIVE EXPANSION 

In this appendix, we carry out a derivative expansion of the effective action 
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S eff = -N Tr In [-D 2 + E) 

= _AT | o °°^Tr( e tD2 -* s ), (Al) 

where 'Tr' represents an integral over space-time and a trace over the internal indices. In the 
second step above, we have used the Schwinger representation. It is an equivalent represen- 
tation at the level of equations of motion and at all orders in the derivative expansion except 
the lowest one. The lowest order yielding the effective potential is handled separately in the 
paper. The Schwinger representation involves an exponential rather than the logarithm and 
is hence better suited for analysis. We thus have to compute the 'trace' 

Tr (e tD2 ' ts ) = J d 4 x-0^e- ikx tr e tD ' z ~ m (e tkx ) 
= J ^(^tr e ^) 2 --(l) 

= / rf 4 x _0_ e -^ tr e -tE +t (2ifcB + ^) (1)> (A2) 

where 'tr' is a trace over the internal indices. One now expands the exponential inside the 
trace and computes different terms to obtain a series representation for the effective action. 
In the literature, to our knowledge, one expands the E term as well. This is not suited for 
our purpose as we intend to keep all orders in E. This suggests that we do a perturbation 
theory in 2ikD + D 2 alone. To this end, we use the following result due to Feynman: 

00 rt ftl 

e -t(H+v) = / dt J dt 2 --- dt n e-^ H (-V)e^-^ H (-V) ■ ■ ■ {-V)e^ H . (A3) 
n=0 Jo Jo 

In our case H = E and —V = 2ikD + D 2 . This leads to an expansion in 2ikD + D 2 . 
Rearranging the terms one obtains a derivative expansion, that is, an expansion in D: 



where a n is of order D 2n . The calculations are quite involved. The result to order D 4 is 



-ts 
1 /■' 



ai = ~2~l ^ e ~ (t_tl)S 2V~' lS , 
a 2 = ~j t f dtM* ~ ii)^e- (t - tl)E F^e- tlS , (A5) 
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where DO for some object O is [-0,0]. The result for ci2 is not complete. However, our 
interest is in its contribution to terms quadratic in with a space-time independent E and 
in this respect it is complete. 

A brief account of the calculations now follows. We write the expansion ( |A3| ) symbolically 

as 

oo 

e -t(H+v) = Y / ^ {-V) ^ (-V) ^ • • - ^ (-V) ^ , (A6) 

n=0 

where again H = S and — V = 2ikD + D 2 . Symbol ^ denotes a 'propagator' of the kind 
exp[— (ti —t2)H]. The beginning and the end of a term of the above kind is indicated by 
black dots. First we note that only terms with an even number of .D's are relevant. Those 
with odd number of Z)'s come with an odd number of k's and their contributions vanish 
after the k— integration. Quite often, we make use of the following reduction to simplify 
results: 

f(t's) ■ ■ ■ Oit^) 1 0(t i+1 ) ■■■= f 1 " drfiU - r) • • • Ofe-O ^ 0{U) • • • , (A7) 

given any function /. In the r.h.s, ti is first absent and we have hence replaced t i+ i by 
ti and so on with t's of higher indices. If, for instance, the function / were absent or is 
independent of ti, this reduction introduces — ti in the r.h.s. Another property we make 
use of to simplify results is the presence of an overall trace and a space-time integral that 
lets us rearrange terms in some expressions. 

Now we come to the calculations. At the lowest order we have a Q = = exp(— tE) 
giving us the effective potential. At the next order, 

a% = 2ikD -o 2ikD oi + ^ D 2 . (A8) 

The first term can be simplified, 

2 

<-^2ikD^2ikD^ = -Ak^r^ D v ^> = — ^D^D^ , (A9) 

t 

where we have replaced k^ky by b^ v j{2t) as the two would yield identical results after 
k— integration. If T> represents the action [D, O] for any O immediately next to it, one 
easily verifies that 
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«-o D r ^ D / ~x-» = «~\-' T> D r ^-» + «~o 1 .D 2 '~\-» 

= *-^V<-^D<-^ + (t-ti)«-x-,D 2 -x-. . (A10) 



Alternately 



*~o D , ~\_> nj = — *-\_, Z}"D '"X-^ 1 '~v_» + <~o Z) 2 o., 1 o» 

= - «-\_, X> 2 1 - ^Oo-Doj + t 1 ^_,_D 2 ^,. (All) 

Adding the two results, 

2*-^D<~^D<-^ = -*-^V 2 +t*-^D 2 <-^ . (A12) 

Putting these together, we have 

t 2 2 v ; 

which is the result quoted earlier. The second step follows from the sum of 

«~o T> 2 '-o 1 u = ti "D 2 and 

«-x_, D 2 -x_, 1 = 1 -x_, £> 2 -x_, = (t-t^^V 2 ^. (A14) 

The next coefficient a>2 can be computed along similar lines. The result quoted earlier 
is obtained by keeping only terms quadratic in with a space-time independent S. This 
simplifies the calculations. The leftmost and rightmost D^s get replaced by — iA^. This 
gives two A^s already so that those D^s in the middle get replaced by <9 M . The contributions 
and their simplified results are 

^D 2 ^D 2 ^ = *^d- A^d- A^* 

2 

*~o D 2 '-o 2ifcD 2i/cD = — (ti-t 2 )^9-A^9-A^ 

2 

«-o 2iA;D ^ D ^ 2ikD = — (t 1 — £ 2 ) *^ c^^L ^ c^A,, r " v -» 

2 

«-\_, 2ifc.D '-o 2i/cZ? <-\_, L> 2 o» = — - t 2 )^ d ■ A^ d ■ A^ 

t 



2i/cD 2i/cD ^ 2i/cD 2i/cD nj = 16k^k u k p k a f^A^ ic^ ^ i<9 p ^ A CT nj . (A15) 
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Here and in the following, a d immediately left to an A as in d ■ A or d ll A u acts only on that 
A, that is, d ■ A = d^A^) for instance. Note that the replacement 

kfJZykpka > (5p U 5p<j + SppSua + ^^cr^up) (A16) 

simplifies the last contribution to 

-(h - t 2 f (2^d-A^d-A^ + ^ dpA u ^ d fJl A v . (A17) 



Adding all the contributions, one gets 



2 



1 - -(ti - t 2 ) 



2 



2 



*^d- A^d- A^-» - -(ti -t 2 ) 



1 



i--(h-t 2 ) 



*-^d ll A v '-^d ll A v '-^. 

(A18) 

Further simplification is possible due to 

f(t 1 -t 2 )^0^0^ = tJit-t^O^O^ (A19) 

and 

f{h-t 2 )^0^0^ = (t-t 1 )f(t 1 )*^0^0= (t-t^fit^O^O^. (A20) 

given any object O. In our case f(t — ti) = f(ti) so that adding and dividing by 2, we get 

f{t x -t 2 )^0^0^ = t -f{t 1 )0^0^. (A21) 

This simplifies the total contribution to 

±-(t-2t 1 ) 2 d- A^d- A^ - -Mt-t^A^dpA^. (A22) 

Note that if one were to work in the Lorentz gauge d ■ A — the first term will vanish and 
the kinetic terms for the gauge fields will arise from the second term. We will regard the 
second term to be a part of the following gauge invariant combination: 

a 2 = -^fidtMt-tJF^e-tt-^F^e-* 1 *. (A23) 

However, the (d ■ A) 2 term that this generates does not agree with what we obtained. 
This is to be expected since there are other gauge invariant combinations, for instance 
X> 2 e~(* _tl ) s D 2 e _ * lE ; , that could generate them. 
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